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The anisotropic two-orbital Hubbard model with different bandwidths and degrees of frustration in 
each orbital is investigated in the framework of both single-site dynamical mean-field theory (DMFT) 
as well as its cluster extension (DCA) for clusters up to four sites combined with a continuous-time 
quantum Monte Carlo algorithm. This model shows a rich phase diagram which includes the appear- 
ance of orbital selective phase transitions, non-Fermi liquid behavior as well as antiferromagnetic 
metallic states. We discuss the advantages and drawbacks of employing the single-site DMFT with 
respect to DCA and the consequences for the physical picture obtained out of these calculations. 
Finally, we argue that such a minimal model may be of relevance to understand the nature of the 
antiferromagnetic metallic state in the iron-pnictide superconductors as well as the origin of the 
small staggered magnetization observed in these systems. 



PACS numbers: 71.10.Fd,71.27.+a,71.30.+h,71.10.Hf 

I. INTRODUCTION 

The multiorbital Hubbard Hamiltonian has been pro- 
posed 1 to be a suitable model for the description of a 
large variety of correlated electron transition metal com- 
pounds including nickelates, vanadates, and more re- 
cently iron pnictide and chalcogenide superconductors 2 . 
While new interesting phases like orbital-selective metal- 
insulator phase transitions, coexistence of paramagnetic 
with antiferromagnetic phases or coexistence of Fermi liq- 
uid with non-Fermi liquid behavior have been predicted 
for this model, there is still a strong debate regarding the 
true nature of the phase diagram. The reason for that 
lies in the complexity of the multiorbital model which 
includes orbital degrees of freedom, crystal field splitting 
effects as well as the Hund's rule exchange coupling and 
also in the nature of the approaches used in order to in- 
vestigate the model. In this work, we shall present a com- 
parative investigation of the properties of the anisotropic 
two-orbital Hubbard model within the single-site dynam- 
ical mean field approximation (DMFT) 3-5 as well as the 
dynamical cluster approximation (DCA) 26 where spatial 
fluctuations are included. 

We will focus on three major phenomena related to the 
multiorbital Hubbard model, (i) Orbital selective phase 
transitions (OSPT) which have been discussed, for in- 
stance, in Ca 2 - a; Sr 2 ;Ru04 6 ' 7 and recently in iron-based 
superconductors 8 ' 9 in the framework of DMFT and den- 
sity functional theory (DFT) calculations. Within single- 
site DMFT it has been shown that the presence of OSPT 
is dependent on various effects like the existence of in- 
equivalent bandwidths at integer filling 10-13 , the presence 
of crystal field splitting at non-integer filling 15 , the ab- 
sence of interorbital hybridization 14 as well as the appear- 
ance of different number of degenerate orbitals 16 . Nev- 
ertheless, it is unclear how these features are modified 
if spatial fluctuations are included in the calculations. 



For instance, the single-site DMFT cannot capture the 
Slater physics caused by nesting effects as is the case 
in the spin-Peierls system TiOCl 17,18 where dimeriza- 
tion is important, (ii) Appearance of non Fermi- liquid 
(FL) behavior. Non-FL behavior is present in the single- 
site DMFT calculations for the multiorbital 20,21 Hubbard 
model and in the four-site cellular-DMFT (CDMFT) 
studies for the one-orbital Hubbard model 22,23 . Based 
on the mentioned calculations, the mechanism of non-FL 
behavior is attributed to the ferromagnetic correlation 
between orbitals (Hund's coupling) in the case of multi- 
orbital systems, while intersite antiferromagnetic corre- 
lations caused by the Coulomb repulsion are suggested as 
the driving mechanism in the one-orbital case. In an at- 
tempt to consider the effect of spatial fluctuations in the 
multiorbital Hubbard model, we investigated recently the 
mechanism of non-FL behavior 24,25 within DCA. This 
study is developed further in the present work, (iii) Fi- 
nally, we investigate the competing effect of frustrated or- 
bitals in a multiorbital Hubbard model. Such a study has 
been motivated by the intensively discussed controversy 
about the origin of the observed reduced magnetic mo- 
ment in Fe-based superconductors 8 and the possible rela- 
tion to the multiorbital nature of the system. Bandstruc- 
ture calculations together with downfolding for various 
families of iron-based SC compounds show that weakly 
frustrated d xy , d z i and highly frustrated d yz / zx , d x i_ y i 
orbitals coexist 32,33 . We will show in what follows that 
competition of frustrated with non-frustrated orbitals in- 
duces an antiferromagnetic metallic phase with a reduced 
ordered magnetic moment 8,25 , as observed in iron-based 
superconductors. 

In section II we present the model and numerical meth- 
ods considered which include the single-site DMFT, DCA 
and the continuous-time QMC algorithm as impurity 
solver. In section III, we discuss the appearance of orbital 
selective phase transitions and non-Fermi liquid behav- 
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ior for two models: (I) the half-filled two-orbital Hubbard 
model with unequal bandwidths on the square lattice in 
the framework of a four-site DCA paramagnetic (PM) so- 
lution and (II) the half-filled two-orbital Hubbard model 
on the Bethe lattice with one frustrated orbital and one 
unfrustrated orbital using the single-site DMFT with AF 
ordering. In this model we find an antiferromagnetic 
metallic phase as well as a reduced staggered magneti- 
zation as a consequence of competing frustration effects. 
Finally, in section IV we summarize our findings. 



II. MODEL AND METHODS 

We consider the anisotropic two-orbital Hubbard 
model with nearest- and next-nearest- hopping matrix 
elements: 

H = — t m cl m(J Cj m v — t m c im<j C j'ma 

(ij)ma- (ij')ma 

+ U n im1 -n im i + ^2(U' - S rja >J z )n ila n l2 a'(l) 

where Ci ma (c im(7 ) is the annihilation (creation) operator 
of an electron with spin a at site i and orbital m. t m (t' m ) 
is the hopping matrix element between site % and nearest- 
neighbor site j (ncxt-nearest-neighbor site /). U and 
U' are intraorbital and interorbital Coulomb repulsion 
integrals, respectively, and J z rii\„nii a is the Ising-type 
Hund's rule coupling term. In our calculations we set 
V = ^ and J z = ^ in order to fulfill the condition 
U = U' + 2J Z and evaluate the model with the single- 
site DMFT and four-site DCA methods combined with a 
continuous-time QMC algorithm. 

Even though DCA can be considered as superior to 
the single-site DMFT method, both approaches are faced 
with certain advantages and drawbacks. While the 
single-site DMFT ignores spatial fluctuations, a notori- 
ous fermionic sign problem present in the QMC calcu- 
lations disappears even in the frustrated systems with 
next nearest-neighbor hopping t' . On the other hand, 
while short-range spatial fluctuations are considered in 
the four-site DCA approach, the QMC calculations en- 
counter a serious fermionic sign problem for the frus- 
trated model 28 . For these reasons, we first consider the 
anisotropic two-orbital model with different bandwidths 
t m and no frustration t' m — on the square lattice at 
half-filling within the four-site DCA method (model (I) 
above). Although a mild fermionic sign problem with av- 
erage sign around (s) = 0.98 is observed, we believe that 
the results are numerically reliable. 

In the DCA method the reduced Brillouin zone (BZ) 
is determined by dividing the total BZ with respect to 
the number of sites and the self-energies for the differ- 
ent cluster momenta are calculated - analogous to the 
DMFT equation - on the reduced BZ by considering 
the same assumption as in DMFT where the self-energy 
is constant in the (reduced) Brillouin zone. The DCA 



self-consistent equations for each cluster momentum with 
K = (0, 0), (0, 7r), (tt, 0) and (w, tt) are given as 

G a (K, ioJn) = ^ T J2 — -Z 1 v rw ■ W 

N ^ iu) n + fi - e K+R - S a (K, ioj n ) 



where N is the number of sites, n is the chemical poten- 
tial, e K+] £ is the dispersion relation, u> n is the Matsub- 

ara frequency, and the summation over K is calculated 
in each reduced Brillouin zone. 

Next, we will study the anisotropic frustrated two- 
orbital Hubbard model on the Bethe lattice at half-filling 
with t m and t' m hopping matrix elements (model (II) 
above) by employing the single-site DMFT with an anti- 
ferromagnetic bath. The Neel state, in which the magne- 
tization of A and B sublattices are in opposite directions, 
is included in the model and the DMFT self-consistency 
equations are given as 

G 0,A,cr = + M — tm^B,a ~ t'm G A,a, (3) 
G 0^B,a = iuJ n + fJ- ~ t 2 m G A,a - t'^G B ,a, (4) 

The sign problem is absent in the single-site impurity cal- 
culations, and the converged impurity Green's functions 
in both cases are determined after several iterations. 

Two possible continuous-time (diagrammatic) QMC 
algorithms for solving the DMFT and DCA equations 
have been proposed. One is based on an expansion 
in the impurity-bath hybridization A 34,35 term - also 
called strong-coupling approach - and the other one is 
based on an expansion in the interaction term (/ 36_38 ) de- 
noted as weak-coupling approach. Unlike a determinant 
QMC method 39 which includes the artificial Trotter error 
caused by discretization of the imaginary time, this error 
is not present in the continuous-time QMC algorithms 
except as numerical noise. In addition, since the matrix 
sizes of both diagrammatic QMC methods are smaller 
than that of the determinant QMC method, it is possi- 
ble to access the low temperature regime. For our study 
we employ the weak-coupling QMC method. The reason 
for this choice lies in the computational effort. While the 
strong-coupling QMC method for the single-site impurity 
problem of multiorbital systems has advantages like the 
possibility of considering the full Hund's rule coupling 
term without any fermionic sign problem and also the 
calculation at very low temperatures in strong-coupling 
regions, the numerical effort of such an approach in- 
creases rapidly with the number of cluster sites N as 2 2N . 
Therefore, in the multiorbital cluster systems the strong- 
coupling QMC method is computationally extremely ex- 
pensive. On the other hand, the weak-coupling QMC 
method contains less numerical burden because the ma- 
trix size increases only linearly with N. 
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III. RESULTS 

A. Orbital selective phase transitions and 
non-Fermi liquid behavior in the four-site DCA 
calculations 

In this section, we present results which demonstrate 
the appearance of OSPT and non-FL behavior in the 
half-filled two-orbital Hubbard model with different hop- 
ping t m in each orbital m on the square lattice (model 
(I)) by employing the four-site DCA calculations with 
paramagnetic bath. The hopping t m of the narrow and 
wide bands are t\ = 0.5 and t 2 = 1-0, respectively, and 
the temperature considered is T/t 2 =0.1. In order to 
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FIG. 1. (a) The quasiparticle weight 2k at K = (ji, 0) as a 
function of U/t 2 for the narrow and wide bands. The spectral 
functions at T/t 2 = 0.1 for (b)-(c) U/t 2 = 2.0, (d)-(e) U/t 2 = 
2.8, and (f)-(g) U/t 2 = 3.4 in the four-site dynamical cluster 
approximation (DCA) calculation. 

investigate the existence of the OSPT, we first explore 
the quasiparticle weight zk as a function of U/t2, which 
is given as zk ~ [l — ImS ^' ItJ °) ] 1 j where luq is the low- 
est Matsubara frequency and K is the momentum sector. 
The effective mass ratio, which is defined as z^ 1 , gives in- 



formation on correlation effects. As is well known, corre- 
lation is the driving force for a transfer of spectral weight 
from the coherent peak at the Fermi level (oj/t 2 = 0) 
to incoherent peaks. This means that the quasiparticle 
weight zk at the Fermi surface K = (tt, 0) should be sup- 
pressed in a Mott insulating state due to enhancement of 
the effective mass. The results are shown in Fig. 1 (a). 
The quasiparticle weight zk at K = (jr, 0) in the narrow 
band is suppressed around U/t 2 — 2.8, while it remains 
large in the wide band. We can denote this feature as 
an orbital-selective phase where metallic and insulating 
states coexist between U/tz = 2.8 and 3.2. In order to 
observe the OSPT more clearly at momentum sectors 
K = (0, 0) and (tt, 0), we plot in Figs. 1 (b)-(g) the spec- 
tral functions for U/t 2 — 2.0, 2.8 and 3.4, where the 
maximum entropy method was used for analytical con- 
tinuation. The spectral functions at K = (0, 0) show a 
band-insulator, while at K = (71", 0) they indicate metal- 
lic states in the weak-coupling region U/t 2 = 2.0 in both 
bands like in results obtained from non-interacting four- 
site DCA in the single-band Hubbard model 40 . Upon 
increasing the interaction we find that the spectral func- 
tions indicating a band insulator at K = (0, 0) are al- 
most unchanged. On the other hand, we observe that for 
U/t 2 = 2.8 the gap is completely opened at K = (tt, 0) in 
the narrow band, while a finite density of states around 
the Fermi level remains in the wide band, clearly show- 
ing the coexistence of an insulator and a metal in this 
two-orbital system as displayed in Figs. 1 (d)-(e). In the 
strong-coupling region with U/t 2 = 3.4 in Figs. 1 (f)- 
(g), both phases at K = (ir,0) are insulating and the 
gap in the narrow band is larger than that in the wide 
band. This means that between U/t 2 = 2.8 and 3.4 a 
Mott metal-insulator transition at K = (tt, 0) occurs in 
the wide band, while the narrow band is already in the 
insulating state; at the same time, a band insulator is 
sustained at K = (0, 0), regardless of the strength of the 
interaction U/t 2 . 

Next, we plot the self-energies as a function of real fre- 
quency in Figs. 2 (a)-(h) to show the FL, non-FL and 
Mott insulating behavior clearly. The analytical contin- 
uation of the self-energy suggested by Wang et al. 41 has 
been used. Here we briefly introduce the method: at high 
frequency oj n the self-energy at half-filling should behave 
as Ti(iuj n ) ~ . In order to perform the Fourier trans- 
formation of Yi(iu> n ), we need to take into account the 
correction of the high frequency tail as 
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where I e"^" T+ ^ = -0.5 and N h is the high- 

est number of Matsubara frequencies. Here, S(r) is in- 
serted into the maximum entropy equations to obtain 
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FIG. 2. Self-energies at the Fermi surface K = (n, 0) as a 
function of real frequency u)/t 2 at T/t,2 = 0.1 for the narrow 
and wide bands in the four-site DCA calculations: (a) and (b) 
U/t 2 = 1.4, (c) and (d) U/t 2 = 2.2, (e) and (f) U/t 2 = 2.6, 
and (g) and (h) U/t 2 = 3.0. 



ImS(o;), where u is the real frequency. ReE(u;) is then 
calculated by the Kramers-Kronig relation. 

In the weak-coupling region, U/t 2 — 1.4, Figs. 2 (a)- 
(b) show that the imaginary part of the self-energy ap- 
proaches zero as lo — > for both orbitals indicating FL 
states. As the interaction is increased, a non-FL state 
in the narrow band coexists with a FL state in the wide 
band as shown in Figs. 2 (c)-(d). By further increasing 
U jt 2 , the narrow band shows a Mott insulating behav- 
ior where ImS(w) diverges at uj — (Fig. 2 (e)), while 
the wide band is still metallic with non-FL behavior (see 
Fig. 2 (f)). In the strong-coupling region U/t 2 = 3.0 
in Figs. 2 (g)-(h), divergences of the imaginary part of 
the self-energy in both orbitals indicate a complete Mott 
insulating state. 



B. OSPT and small magnetic moment in the 
single-site DMFT calculation 

We proceed now to the analysis of the anisotropic two- 
orbital Hubbard model with frustrated and unfrustrated 
orbitals at half-filling using the two-sublattice single-site 



DMFT with antiferromagnetic bath. In our calculations, 
the temperature is set to T/t = 1/16. The next near- 
est neighbour hopping is t'/t = for the unfrustrated 
orbital and t'/t = 0.65 for the frustrated orbital with 
t = 1.0. We first evaluate the OSPT behavior of this 
model with the single-site DMFT, as we did within the 
four-site DCA approach in the former section. We ob- 
serve that a band insulator is present since the source 
of the gap opening is not only the Coulomb interaction 
but also the magnetization of the AF solution. There- 
fore, the change of the effective mass - which is related 
to the quasiparticle weight z - is very small in contrast 
to the case of a Mott insulator. Therefore, we employ 
the DOS at Fermi level P (uj = 0) to identify the OSPT 
as is shown in Fig. 3. The results were obtained through 
analytical continuation applying the maximum entropy 
method and we performed around 10 7 QMC samplings 
and checked the results through a Pade approximation 
in order to ensure the reliability of the results. The gap 
opening is clearly observed around U Cl /t = 1.9 in the un- 
frustrated orbital, while the metallic state is still present 
in the frustrated orbital. As the interaction is increased, 
a metal-insulator transition in the frustrated orbital ap- 
pears around U C2 /t = 2.9. Between U Cl /t < U/t < U C2 /t 
we find a coexistence region of metallic and insulating 
states as a clear evidence of OSPT behavior. Compar- 
ing our results of Fig. 3 with those obtained with the 
single-site DMFT with a paramagnetic solution in the 
two-orbital Hubbard model with different bandwidths 10 , 
we find that in both cases the spectral densities p(ui = 0) 
decay suddenly to zero in the narrow (unfrustrated) band 
and the values in the wide (frustrated) bands are similar 
around U Cl . The difference between these two results is 
that p(cj = 0) for the frustrated orbital of our system 
also drops to zero at U C2 , while it decreases smoothly 
in the single-site DMFT calculation with paramagnetic 
solution . We believe that this is due to the different 
mechanisms of the gap opening in each system. 

The present model is also relevant for understanding 
the magnetic behavior of the iron-pnictide superconduc- 
tors. Since the discovery of high-T c superconductivity in 
iron pnictides in 2008 2 , a lot of studies have been per- 
formed to find the mechanism of pairing in these mate- 
rials 42 ~ 44 . Unlike the high-T c cuprate superconductors, 
whose parent compounds are antiferromagnetic insula- 
tors, an ordered antiferromagnetic metallic state with 
small magnetic moment is detected in iron pnictides. 
Density functional theory calculations on iron pnictides 
overestimate the values of the ordered magnetic moment 
compared with experimental results 45-48 . Only by artifi- 
cially considering a negative U or by employing a specific 
double counting correction within LDA+[/ 49 ' 50 can these 
values be reduced . Even though some DFT+DMFT cal- 
culations including quantum fluctuations have been per- 
formed, the mechanism of the small magnetic moment 
is still controversial 5 1-63 . Here, within model (II) we 
discuss a possible mechanism for the small magnetic mo- 
ment in iron pnictides and the presence of a metallic an- 
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FIG. 3. The density of states p(uj — 0) at the Fermi level 
as a function of U /t for the unfrustrated and frustrated or- 
bitals in the single-site dynamical mean field theory (DMFT) 
calculations. 



tiferromagnetic state. 
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FIG. 4. Magnetization as a function of U/t with different 
values of Jz in the single-site DMFT calculations. 

In our former studies, we proposed the coupling be- 
tween frustrated and unfrustrated orbitals in a multi- 
orbital Hubbard model as the possible mechanism for 
the small magnetic moment observed in iron pnictides 8 . 
Based on the minimal two-orbital Hubbard model pre- 
sented above, we found 8 a continuous metal to insulator 
phase transition as well as a continuous magnetic phase 
transition for an Ising-Hund's rule coupling of J z = U /A. 
The reduced ordered magnetic moment is shown in Fig. 4 
as a function of U/t. We observe at a critical U/t a 
smooth increase of the magnetic moment, which is in 
contrast to the first order phase transition behavior ob- 
served in the one-orbital frustrated Hubbard model 64 . 
In order to show that the results are independent of the 
choice of J z we also show in Fig. 4 the magnetization 
for J z = 7[//24, U/8 and [7/4 with the constraint of 
U = U' + 2J Z . As the interaction is increased, the mag- 



netization in all cases increases smoothly and in the in- 
termediate interaction regime small magnetizations are 
observed. 



IV. SUMMARY 

In the present work we have explored the properties of 
two models. (I) The two-orbital Hubbard model with 
unequal bandwidths on the square lattice by employ- 
ing the four-site DCA and (II) the two-orbital Hubbard 
model with a frustrated and an unfrustrated orbital on 
the Bethe lattice by using the single-site DMFT method. 
The DCA and DMFT equations were solved using a 
continuous-time QMC algorithm. The above choice of 
systems and methods explicitly avoids the fermionic sign 
problem in the QMC calculations. We find that both 
models show a rich phase diagram as a function of inter- 
action strength which includes the appearance of orbital 
selective phase transitions, non- Fermi liquid behavior as 
well as antiferromagnetic metallic states. 

Using the four-site DCA method with a PM solution 
for the two-orbital Hubbard model with different band- 
widths (model (I)) we find that while the system behaves 
as a band insulator in the K=(0,0) sector of the Brillouin 
zone, it undergoes an orbital selective phase transition in 
the K=(7r,0) sector. The system in the K=(7r,0) sector 
shows FL behavior for both orbitals in the weak-coupling 
region. As the interaction is increased, a non-FL state for 
the narrow band coexists with a FL behavior of the wide 
band. In the intermediate-coupling regime the non-FL 
behavior is changed into a Mott insulating state where 
the electrons are completely localized in the narrow band, 
while the metallic state with non-FL behavior is still ob- 
served in the wide band. In the strong-coupling region 
the spectral weight in each orbital shows a gap at the 
Fermi level, indicating a Mott insulating state in both 
orbitals. The observation of an OSPT in the intermedi- 
ate region is due to the different electron localization in 
the orbitals. The non-FL behavior is also present due to 
weakly ordered states with AF correlations between sites 
and FM correlations between orbitals in the metallic re- 
gions. 

Next, using the single-site DMFT with an AF solu- 
tion, we studied model (II) and found that the critical 
interaction U for a metal-insulator phase transition are 
U Cl /t =1.9 and U C2 /t — 2.9 for the unfrustrated and 
frustrated orbitals, respectively. This is a clear evidence 
of coexistence of two phases between the two critical U 
values. An OSPT is observed in this system since the 
strength of electron localization in both orbitals is dif- 
ferent due to the unequal orbital frustration. The AF 
states are a key to the gap opening and the nature of 
the insulating state is that of a band insulator. Finally, 
such a minimal model may be of relevance to understand 
the nature of the antiferromagnetic metallic state in the 
iron-pnictide superconductors as well as the origin of the 
small staggered magnetization observed in these systems 
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since we measure the magnetization in each orbital and 
find that the magnetization only increases smoothly as a 
function of the interaction strength in the antiferromag- 
netic metallic state. 
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